Observation of the nuclear magnetic octupole moment of 173 Yb from precise 
measurements of hyperfine structure in the 3 P 2 state 
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We measure hyperfine structure in the metastable 3 P2 state of 173 Yb and observe the nuclear 
magnetic octupole moment. We populate the state using dipole-allowed transitions through the 3 Pi 
and 3 Si states. We measure frequencies of hyperfine transitions of the 3 P2 — > 3 Si line at 770 nm 
using a Rb-stabilized ring cavity resonator with a precision of 200 kHz. Second-order corrections 
due to perturbations from the nearby 3 Pi, 3 Po, and 1 Pi states are below 30 kHz. We obtain 
the hyperfine coefficients as: A = -742.11(2) MHz, B = 1339.4(2) MHz, which represent two 
orders-of-magnitude improvement in precision, and C = 0.54(2) MHz. 

PACS numbers: 32.10.Dk,32.10.Fn,32.80.Xx,31.30.Gs 



Observation of the nuclear magnetic octupole moment 
and its influence on the hyperfine structure of atoms 
has remained largely unexplored because of its weaker 
effect compared to the leading magnetic dipole and elec- 
tric quadrupole moments. The most recent measure- 
ments have been on the short-lived 6p 3 / 2 state of the 
one-electron atom 133 Cs [H, 0], but the sub-kHz value re- 
quired similar precision in measuring the intervals. By 
contrast, the long-lived 3 P 2 state in two-electron atoms 
such as Yb and several alkaline-earth metals, with its 
large angular momentum, is a more sensitive probe to 
observe this moment Q. Here, we report precise mea- 
surement of hyperfine structure in the 3 P2 state of 173 Yb, 
and first observation of the magnetic octupole moment 
in this class of atoms. Precise measurement of hyperfine 
structure is also motivated by the fact that its compari- 
son to theoretical calculations plays an important role in 
validating the atomic wavefunctions used in the calcula- 
tions. In this regard, Yb is an important atom because of 
its proposed use in search for a permanent electric dipole 
moment (EDM) 0,i], where comparison to calculation 
@> 13 forms a vital tool in searching for new physics be- 
yond the Standard Model. The 3 P2 state in Yb has po- 
tential applications in more sensitive EDM searches Q, 
and ultra-sensitive magnetometry [9]. 

Measurements on upper levels is an experimental chal- 
lenge because these levels are not directly populated. We 
have earlier solved this problem by using dipole-allowed 
transitions to pump atoms into the metastable 3 P 2 state 
of Yb [8j . In this work, we use the same method to pop- 
ulate this state, and then measure the absolute frequen- 
cies of various hyperfine transitions on the 3 P2 3 <Si 
line. We measure the frequencies with our well-developed 
technique of using a Rb-stabilized ring-cavity resonator 
p}]Hl~3j . We obtain the hyperfine structure constants A 
(magnetic dipole) and B (electric quadrupole) with two 
orders-of-magnitude better precision than previous val- 
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FIG. 1: (Color online) Spectroscopy chamber for optically 
pumping atoms into the 3 P2 state and measuring the 3 P2 — > 
3 Si transition at 770 nm. 



ues, and a 4% measurement of the magnetic octupole 
constant C. We take into account second-order correc- 
tions due to perturbations from the nearby 3 Pi, 3 Po, and 
1 Pi states, and find that they make negligible contribu- 
tion. 

The ground state in Yb is 1 So and therefore has no 
hyperfine structure. But the upper states with J 
have hyperfine levels for the odd isotope 173 Yb, deter- 
mined by the nuclear spin 7 = 5/2 8]. The various tran- 
sitions are accessed in the spectroscopy chamber shown 
in Fig. [TJ It consists of a vacuum chamber with sev- 
eral optical access points maintained at a pressure below 
1CP 8 mbar with an ion pump. The Yb atomic beam is 
generated by resistive heating (to about 400°C) of an un- 
cnriched source. The different laser beams are sent per- 
pendicular to the atomic beam at different points, and 
the green fluorescence (at 556 nm) is collected by two 
photomultiplier tubes (PMT, Hamamatsu R928). The 
556 nm beam, driving the So — > 3 Pi transition, is pro- 
duced by doubling the output of a fiber laser operating 
at 1111 nm (Koheras Boostik Y10). The output power 
of the fiber laser is 0.5 W with a linewidth of 70 kHz. It 
is doubled in an external cavity doubler to give a total 
power of 65 mW. Part of this beam is split into two and 
sent across the atomic beam in counter-propagating di- 
rections to minimize systematic Doppler shifts. The laser 
is locked to the peak center by frequency modulation at 
20 kHz and lock-in detection to generate the error signal. 

The lasers at 680 nm (driving the 3 Pi — > 3 Si transi- 
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tion) and 770 nm (driving the 3 P2 —> 3 Si transition) are 
home-built diode laser systems [lj]. They are frequency 
stabilized using grating feedback to give linewidths of 
order 1 MHz. The 680 beam counterpropagates with 
the locked 556 beam a few cm downstream from the 
556 fluorescence point. It optically pumps atoms into the 
metastable 3 Pq and 3 P 2 states. This reduces population 
in the ground state and the green fluorescence spectrum 
when scanning the 680 beam shows a negative peak. The 
injection current into this laser is frequency modulated 
(at 15 kHz) and it is locked to the peak center. The 770 
beam, which is a further 2 mm downstream and overlaps 
with the same 556 beam, causes the fluorescence to re- 
cover as it repumps atoms from the 3 P2 state back into 
the ground state, from where they can reabsorb the 556 
beam. Thus, the 770 beam measures the population in the 
P 2 state and the fluorescence spectrum shows a positive 
peak. The measurements are done with this laser locked 
to a particular hyperfine transition. For this, the laser is 
current modulated at a slightly different frequency from 
that of the 680 laser, and the same PMT signal is de- 
modulated at this frequency to generate the error signal. 

Representative fluorescence spectra in 174 Yb obtained 
by scanning each of the three lasers have been shown in 
our previous work |8j . Here, we show a typical spectrum 
for a particular hyperfine transition in 173 Yb, but only 
for the 770 laser (Fig. [5]). The solid line is a curve fit 
to a Voigt profile, with the fit residuals shown on top. 
The fit yields a linewidth of about 36 MHz. For com- 
parison, the natural linewidth of the 3 P2 —> 3 Si line 
is 12 MHz. Considering that the fluorescence spectrum 
measures the final population in the ground state, this in- 
crease in linewidth is reasonable. From the number of hy- 
perfine levels and the electric dipole selection rules, there 
are 10 possible transitions between these states. The 
spectra for the different transitions are similar to the one 
shown, but appear with different signal-to-noise ratios 
(SNR) depending on the individual transition strengths. 
The SNR (of about 15 from the fit residuals) for the spec- 
trum shown is neither the best nor the worst, but in all 
cases it is good enough to get a strong error signal to lock 
the laser. 

Our frequency-measurement technique using the Rb- 
stabilized ring-cavity resonator has been described ex- 
tensively in earlier work [ToL [l2j , and is reviewed here for 
completeness. It relies on the fact that the frequency of 
the 780-nm D 2 line in 87 Rb (5Si/ 2 o 5P 3 / 2 transition) 
is known with an accuracy of 6 kHz [15j. A diode laser 
locked to this line is used as a frequency reference. A 
second laser is in turn locked to the unknown transition. 
The two lasers are coupled into an evacuated ring-cavity 
resonator. An acousto-optic modulator (AOM) placed 
in the path of one of the two lasers is used to produce 
a small frequency offset (of order 100 MHz) so that the 
cavity is in simultaneous resonance with both laser fre- 
quencies. Thus the ratio of the unknown frequency to 
the reference frequency is just a ratio of two integers 
(i.e. the respective cavity mode numbers) combined with 
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FIG. 2: (Color online) Green fluorescence signal as the 770 
laser driving the 3 P2 3 Si transition is scanned. The solid 
curve is a fit to a Voigt profile with the fit residuals on top. 
The spectrum shown is for the F = 7/2 — > 5/2 hyperfine 
transition. Other transitions give similar spectra but with 
varying signal-to-noise ratio. 



the AOM offset. The procedure to determine the unique 
mode-number combination for the two lasers has been 
described earlier [lj|]. We have used this technique to 
measure frequencies of transitions in the range of 670 nm 
to 895 nm, i.e. ±100 nm away from the reference wave- 
length. 

Our technique is particularly well-suited to the mea- 
surement of frequency differences (as used in the mea- 
surement of hyperfine intervals). This is because several 
sources of systematic errors cancel when taking the dif- 
ference. Systematic errors can be classified under three 
categories: errors related to locking the reference laser, 
errors related to locking the unknown laser on the Yb 
line, and errors due to the measurement cavity. Errors 
related to the reference laser are the same for all measure- 
ments and therefore cancel in the difference. Similarly, 
errors due to dispersion inside the cavity and at the cavity 
(multilayer dielectric) mirrors do not affect the difference. 
Effects of geometric misalignment into the cavity result in 
the excitation of higher-order modes, which can pull the 
lock point of the fundamental mode. To the extent that 
the misalignment is the same for all measurements, this 
will again not affect the difference. Thus all the sources 
of systematic error in the difference frequency are related 
to how well we can lock the 770-nm laser to a particular 
Yb transition and how well we can lock the cavity to the 
laser, as discussed below. 

Collisional shifts of the transition frequency are min- 
imal because we use an atomic beam. A systematic 
Doppler shift of the peak center will occur if the angle be- 
tween the atomic beam and the laser beam is not exactly 
90°. However, this will be the same for all transitions 
and will not shift the interval. Another potential source 
of systematic shift in the peak position is due to line- 



3 



shape asymmetry that might occur because of selective 
driving into Zeeman sublevels (shifted in the presence of 
stray magnetic fields) or radiation pressure. The first ef- 
fect is minimized by using linearly polarized light so that 
the Zeeman sublevels are excited equally about line cen- 
ter. The maximum Zeeman shift among all transitions 
in the presence of a stray field of 10 mG is about 70 kHz. 
By studying the fluorescence lineshape in Fig. [3J we con- 
clude that this error is smaller than 70 kHz. Thus, the 
main source of error in our measurement is determined 
by how well we lock to peak center (split the line) given 
the linewidth of about 40 MHz and the SNR. In our re- 
cent work on measuring isotope shifts in the 556-nm line 
of Yb [13], the linewidth was about 6 MHz and the lock- 
ing error was 30 kHz. With the 6 times larger linewidth 
here (Fig. [5]), we estimate the error to be 200 kHz (or 
equal to splitting the line by about 1 part in 200), which 
is larger than all other sources of error. As we will see be- 
low, we have a good experimental handle on this estimate 
by measuring the same transition with two different lock 
points of the reference laser, one from the F = 1 — y F' 
hyperfine manifold and the other from the F = 2 — > F' 
manifold. 

We have measured the frequencies of various hyper- 
fine transitions by locking the 770-nm laser to different 
peaks. For each measurement, the time constant in the 
frequency counter (with a timebase stability of 10 -8 ) was 
set to 10 s. Then a set of 40 independent measurements 
was made and the average determined. For each transi- 
tion, the frequency was measured with the reference laser 
on either a F = 1 — > F' transition or a F — 2 — > F' tran- 
sition. The resulting change in the frequency is on the 
order of 6.5 GHz, and is known with < 10 kHz precision 
from hyperfine measurements in 87 Rb (l5j . 

The measured frequencies are shown in Table HI The 
error in each value is 0.2 MHz. In order to highlight the 
above fact that many sources of error cancel in the dif- 
ference, we show the frequencies as offset from the first 
value. We have also verified that if we make an error in 
the mode number and change it by ±1, the absolute fre- 
quency of all transitions changes by about 12 MHz but 
there is only negligible change in the offsets listed in the 
table. Our error estimate is reasonable because the two 
values for each transition, which implies a different set 
of cavity mode numbers and complete re-optimization of 
all the feedback loops, overlap quite well. In addition, 
the maximum standard deviation in each set of 40 mea- 
surements (from which the average is determined) is only 
170 kHz. 

In order to fit the measured intervals to the hyperfine 
coefficients, we first need to know the energy shift of an F 
level due to the hyperfine interaction. From perturbation 
theory, the shift is given by 

w F = wg> + wP , (1) 

where the first-order shift has the progressively weaker 
magnetic dipole (A), electric quadrupole (B), and mag- 
netic octupole (C) terms. The prefactors multiplying 



these coefficients for each F level in the 3 P2 state are 
listed in Table M 

The second-order shift in Eq. ([1} arises mainly be- 
cause of perturbations from the nearby 3 Pi, 3 Pq, and 
1 Pi states [16j . The shift is a function of the parameters 
a s, »3/2j ^3/2; £) an d T). Here, a's and b's are defined 
in terms of the hyperfine coefficients of the respective 
states, and £ and r\ are defined in terms of the ratios 
of the off-diagonal to diagonal matrix elements for the 
dipole and quadrupole operators respectively. Using pre- 
viously measured values for the hyperfine constants in 
Refs. jli"| and [l3|, and the interaction terms £ = 1.02 
and rj = 1.26 from Ref. 17], we can calculate the second- 
order corrections listed in Table [TT1 

The procedure to extract the coefficients is now 
straightforward. The shift of each level in terms of A, 
B, C, and the second-order correction is known. This 
gives a set of 6 equations corresponding to the 6 mea- 
sured transition frequencies. We then do a least-squares 
fit to these equations with the hyperfine coefficients as 
fit parameters. In order to see the relative importance 
of C and the second-order correction, we have calculated 
the coefficients with and without these parameters. The 
results are shown in Table Hill The difference from not in- 
cluding the second-order corrections is negligible, which 
is to be expected since the corrections are much smaller 
than our overall error bars. However, not including C 
makes a large difference: of 15c (combined) in A and of 
Her (combined) in B. This shows that there is compelling 
evidence for C from our measurement, and the extracted 
value of 0.54(2) MHz demonstrates that its value lies be- 
tween 0.46 and 0.60 MHz with 95% confidence. We have 
also tried to see if there is a signature of the nuclear 
electric hexadecapole coefficient D, and find a value of 
0.3(30) kHz consistent with zero. 

Our values are compared to previous values in Table 
IIV1 The most recent experimental value is a 1992 mea- 
surement by Maier et al. [H[ . The two sets of values are 
consistent but our error bars are about 100 times smaller. 
An earlier 1962 measurement by Ross and Murakawa [l9| 
gives values in cm -1 with no error bars. Their values 
appear close but we cannot judge the overlap without er- 
rors. Not surprisingly, none of these values are sensitive 
to the octupole coefficient C. There is also a theoreti- 
cal calculation in 1999 by Porsev et al. [2(|, which yields 
values within a couple of MHz (less than 0.5%) of our 
values. 

In conclusion, we have made a precision measurement 
of hyperfine structure in the 3 P2 state of 173 Yb, and see 
an unmistakable signature of the magnetic octupole con- 
stant C. The frequencies of the 3 P2 — > 3 Si transition at 
770 nm are measured using a Rb-stabilized ring-cavity 
resonator with an accuracy of 200 kHz. Second-order 
corrections are negligible at this level of precision, and 
do not affect the value of C . We plan to complete sim- 
ilar measurements in the other odd isotope, 1 Yb, and 
thereby get a handle on higher-order effects such as the 
Bohr-Weisskopf effect and the hyperfine anomaly. 
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TABLE I: Measured frequencies for various hyperfine transitions of the 3 P2 —5- 3 Si line with the reference laser locked to either 
an F = 1 — } F' or an F = 2 — Y F' transition. The values are given as offset from the first transition. The numbers in brackets 
are la errors in the last digit. 



Ref. Laser 






Frequency (MHz) 






lock point 


5/2 ->• 3/2 


5/2 -> 5/2 


5/2^7/2 7/2^5/2 


7/2 7/2 


7/2 9/2 


1 ->■ F' 

2 ->• F' 




0.21(20) 


2523.43(20) 
2523.56(20) 


5360.52(20) -4075.69(20) 
5360.60(20) -4075.65(20) 


-1238.38(20) 
-1238.46(20) 


1193.92(20) 
1193.97(20) 



TABLE II: Prefactors multiplying the hyperfine coefficients of 
the first-order shift for each F level in the 3 P2 state of 173 Yb. 
The last row is the second-order correction to the shift in 
MHz. 



F 


1/2 


3/2 


5/2 


7/2 


9/2 




-7 


-11/2 


-3 


1/2 


5 




7/10 


1/4 


-1/4 


-17/40 


1/4 




-42/5 


12/5 


27/5 


-22/5 


1 







-0.028 


-0.031 


0.049 
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TABLE III: Calculated values of the hyperfine constants for 
the 3 P2 state in 173 Yb. The second row is calculated without 
including the second-order correction and the third row is 
calculated without the magnetic octupole constant C. All 
values in MHz. 





A 


B 


C 


Corrected 


-742.11(2) 


1339.4(2) 


0.54(2) 


Uncorrected 


-742.11(2) 


1339.4(2) 


0.55(2) 


Without C 


-742.53(2) 


1342.5(2) 





TABLE IV: Comparison of hyperfine constants for the 3 P? 
state in 173 Yb from this work to previous values. 



A 


B 


C 


Reference 


-742.11(5) 


1339.4(5) 


0.55(4) 


This Work 


-742(5) 


1342(38) 




[18] 


-738 


1310 




[19]° 


-745 


1335 




Theory [20] 



"1962 measurement with values given in cm 1 and no error bars. 
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